Robust d-wave Pairing Correlations in the Heisenberg Kondo Lattice Model 
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The Kondo lattice model enlarged by an antiferromagnetic coupling Jaf between the localized 
spins is here investigated using computational techniques. Our results suggest the existence of 
a d-wave superconducting phase close to half-filling mediated by antiferromagnetic fluctuations. 
This establishes a closer connection between theory and heavy fermion experiments than currently 
provided by the standard Kondo lattice model with Jaf=0. 
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Introduction. Heavy fermions (HF) continue attract- 
ing the interest of the condensed matter community []]]. 
These materials are phenomcnologically described by Do- 
niach's scenario, where localized spins interact with con- 
duction electrons via an exchance interaction J Q. At 
small J, the indirect Ruderman-Kittel-Kasuya-Yosida 
(RKKY) mechanism is expected to induce an antiferro- 
magnetic (AF) state, whereas for large J a paramagnetic 
spin-liquid state emerges [2| . The discovery of supercon- 
ductivity in CeCu2Si2 [H, and subsequently in several 
other HF compounds, unveiled the rich variety of phe- 
nomena present in these strongly correlated electronic 
systems. Currently, HF superconductors are widely con- 
sidered as examples of superconductivity not mediated 
by lattice vibrations 0, Q • 

Several formal theories of HF materials start with the 
Kondo lattice model (KLM) [f|. The development of 
powerful many-body numerical techniques, and the con- 
tinuous growth in computer power, have allowed for non 
mean-field and free of non-controlled parameters inves- 
tigations of the KLM, at least in low dimensional sys- 
tems. These investigations have revealed two potentially 
important problems in establishing a connection between 
KLM and HF phenomenology: (1) Including hole carriers 
away from half-filling, and at large J, the KLM leads to 
a robust ferromagnetic state that is not obviously con- 
nected with states found experimentally Q; (2) more 
importantly, there is evidence that the standard one- 
dimensional (ID) KLM does not present SC tendencies 
close to half-filling [8( . Although these anomalies may be 
caused by the low dimensionality of the Kondo lattices 
studied, they still raise doubts about the full validity of 
the simple KLM to describe heavy fermions. Moreover, 
it would be advantageous for theoretical investigations 
to identify a simple HF model with clear SC tendencies, 
even in low dimensions. The growing number of SC HF 
compounds clearly requires a model paradigm involving 
not only AF and spin-gapped phases, as in the past, but 
including a SC phase as well. 



To address these concerns, and better capture the 
physics of HF systems, we must move beyond the stan- 
dard KLM. In this letter, using computational techniques 
it is shown that the addition of a direct AF coupling Jaf 
between the localized spins considerably alleviates the 
problems mentioned above. Previous investigations had 
already suggested that for UPd2Al3 the AF coupling be- 
tween localized spins is important to understand the SC 
phase 0. Moreover, other authors already remarked the 
importance of Jaf when focusing on the magnetism of 
the HF systems [13, 11, [Hj]. However, our present ef- 
fort goes beyond these previous studies by revealing the 
formation of d-wave symmetric hole-pairs and the devel- 
opment of robust SC correlations when Jaf is incorpo- 
rated to the two-dimensional KLM. This establishes a 
better connection theory-experiment in the HF context 
than provided by the plain KLM. 

Model/Methods. The antiferromagnetic Heisenberg 
Kondo Lattice Model (HKLM) is given by 
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where J>0 is the Kondo coupling between the conduc- 
tion electrons and the local moments, Jaf is the antifer- 
romagnetic interaction between the localized spins-1/2, 
and the hopping amplitude was set to unity to fix the 
energy scale. The rest of the notation is standard. The 
HKLM on a NxL cluster was here investigated u sing 



14] 



the Lanczos technique [13|, and the DMRG method 
under open boundary conditions (OBC) [lil. Il6j|. 

Note that the same HKLM also describes the man- 



ganites if J<0 and the Sj's are assumed classical [17 1 



Investigations of manganitc models have shown that Jaf 
is crucial for the numerical stabilization of experimen- 
tally known phases that otherwise become unstable due 
to the strong ferromagnetic tendencies 17, 18} . Thus, our 
effort also provides a unifying view of HF and manganite 
research regarding the relevance of the Jaf coupling. 
Binding Energies. Our investigations mainly focused 
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FIG. 1: (Color online) The region with binding of holes in the 
J- Jap plane, for different cluster shapes with OBC. The solid 
line corresponds to Ab=— 0.1. 



on pairing and SC correlations. To observe indications 
of pairing in the HKLM ground state, we measure the 
pair-binding energy defined as Ab=E(Q)+E(2)-2E(1), 
where E{1) is the ground state energy in the subspace 
with (NL — I) conduction electrons. If the holes do not 
form a bound state, in a finite system the binding energy 
is positive Ab>0 [131 ]. while in the thermodynamic limit 
it vanishes. On the other hand, if holes do form a bound 
state then Ab<0, and this is indicative that attractive 
effective forces are present, as widely discussed before in 
the context of high-T c investigations [HI . 

In Fig. 1, we show the region in the J-Jaf plane where 
Ab is less than -0.1, for the cluster sizes 1x10, 2x6 and 
for a square lattice 3x3 with OBC. We choose to present 
the region where Ab<— 0.1, as opposed to Ab=0, since 
previous experience in the cuprate context [l3| suggests 
that this procedure takes better into account the size ef- 
fects. The limitation of using small clusters originates in 
the huge Hilbert space of the HKLM, with 8 states per 
site, and in the need to calculate hundreds of points to ex- 
tract comprehensive phase diagrams. However, for some 
selected couplings much larger clusters were considered 
and a thermodynamic limit extrapolation was carried out 
(see Fig. 2 and discussion below). Note that the bind- 
ing regions in Fig. 1 are qualitatively similar, although 
the clusters have different shapes. Thus, the formation 
of hole pairs is a robust effect, suggesting that pair for- 
mation might also appear in the three dimensional phase 
diagram (currently unreachable numerically). Moreover, 
the tendency toward increasing pairing strength found in 
Fig. 1 moving from chains to ladders/square clusters is 
similar to tendencies found in the t-J model for cuprates 
fl3j |. It is interestin g, al so to note that a previous study 
by Sikkeman et al. have found that Jaf induces a 
spin gap phase with no charge gap close to half-filling 
for the one-dimensional HKLM, in agreement with the 
phase diagram present in Fig. 1. Finally, note that val- 
ues of Jaf are expected to be small due to the small 
overlap between nearest-neighbor 4f/5f orbitals, and in 
Fig. 1(c) it appears that hole binding can occur for Jaf 
as small as 0.2. Only a careful two-dimensional size scal- 
ing, beyond our capabilities at present, can answer how 




FIG. 2: (Color online) Pair-binding energies Ab vs. 1/L for 
a chain (a), a two- leg ladder (b), and square clusters (c). (d) 
Ab vs. 1/m for the cluster 6x6, at couplings indicated. 



small Jaf can be to still induce hole pairing. 

In Figs. 2(a)-(c), Ab vs. 1/L is shown for the chain, 
two-leg ladder and square clusters, respectively. There 
are clearly two distinct behaviors. For the couplings that 
present hole binding for the small clusters presented in 
Figs. 1, Ab tends to negative values when L—>oo, while 
for the other parameters Ab^O. Regarding the 6x6 clus- 
ter, even working with up to to=4000 states the energies 
did not converge well. However, as presented in Fig. 
2(d), Ab tends to a negative value for the intermediate 
coupling range, similar as the results obtained for the 
other cluster shapes. All these results indicate that close 
to half-filling, pairing tendencies exist for intermediate 
values of J and Jaf in the bulk of the HKLM. It is also 
interesting to observe the behavior of Ab vs. J for a fixed 
Jaf (or fixing J and varying Jaf)- As present in Fig. 
3(a), — Ab reaches a maximum value at intermediate J, 
important detail to compare with experiments, as dis- 
cussed later. We found that the hole binding is robust in 
the same region where the nearest-neighbor spin-spin cor- 
relations are robust, as shown in the inset of Fig. 3(a). 
There, we also plot the nearest-neighbor hole-hole cor- 
relation [ijj]. The fact that spin and hole correlations 
behave similarly supports our claim that antiferromag- 
netism and hole binding are related. At very large J, 
both quantities are suppressed together. This suggest 
that the origin of the hole-pair attraction is connected 
with AF fluctuations. We believe mechanisms similar to 
those identified in cuprates 



13j, such as hole attraction 



to minimize the "damage" to the AF background, are in 
action in the HKLM as well 

Pair Symmetry. Our exact diagonalization results for 
the clusters 2x2 and 3x3 show that the ground state 
with (2) holes has s-wave (J-wave) symmetry under ro- 
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FIG. 3: (Color online) (a) —As vs. J for the clusters 1x10, 
2x6, and 3x3 with OBC. Inset: The nearest-neighbor spin- 
spin (S(l)) and hole-hole (JV(1)) correlations, measured in 
the center of the 2x6 cluster with 2-holes and Jaf=0.75, as 
function of J. (b) The pair-pair correlation function C(l) vs. 

1 for a system size 2x32 and Jaf=0.75 with 8 holes, for some 
values of J. The straight lines are data fits with the powers 
indicated by the arrows, (c) P vs. J for the cluster 2 x 32 with 

2 and 8 holes and some values of Jap- (d) The spin structure 
factor S(n, n) vs. Jaf for the clusters 1x20, 2x10 and 4x4, 
at J=2 and half-filling. 



tations. Thus, these ground states are connected via a 
pair-creation operator with d-wave symmetr y, s imilarly 
as observed in the Hubbard and t-J models [13J. Then, 
the HKLM predicts <i-wave superconductivity, compati- 
ble with previous results for HF systems [5| . 

Pairing Correlations. To confirm that hole-pair ten- 
dencies are concomitant with robust superconducting 
correlations, for the two-leg ladder we measured the rung- 
rung pair correlation functions C(l) — (^Ai + i where 

the pair operator is defined as A- J = |Q2,j.-Cq |C;2,t, an d 
cix,<t annihilates a conduction electron on rung I and leg 
A = 1,2 with spins a =f,| 0- In Fig. 3(b), C(l) is 
shown vs. distance I, for a fixed Jaf and some values 
of J, using a cluster 2x32 with 8 holes. Similar results 
are obtained with 2 holes. In order to obtain the slope of 
C(l) we fit our data with the function a/L n (n an inte- 
ger). Clearly from this figure, the pair-pair correlations 
C(l) are enhanced at large distances for intermediate cou- 
pling values, i. e. J=3 and Jaf=0.75, with a robust 
power-law decay ~1/Z. It is interesting to note that the 
well known two-le g t- J model also has the same slope 
close to half-filling [2 lL I22I ] . The origin of this effect may 
be due to the fact that the HKLM can be mapped, in a 
limiting case, into the t-J model [ll|. Similar SC tenden- 
cies are also observed in Fig. 3(c), where f=X^i=gC(0 
is shown vs. J, at a fixed Jaf- P is enhanced only for 
non-zero values of Jaf, indicating that superconductivity 



appears only when this interaction is active. Although 
there is no true long-range order in quasi-lD systems, 
our results provide strong evidence that superconductiv- 
ity dominates at intermediate J and Jaf in the HKLM 
close to half-filling. Similar conclusions were reached for 



the t-J model on chains and ladders 

Magnetic Properties. We also investigated the effect of 
Jaf in the magnetic properties of the HKLM. At small 
values of J, due to the the RKKY interaction, antiferro- 
magnetic long-range order is expected, whereas for large 
J a paramagnetic state must emerge. The competition 
between these two states leads to a quantum critical point 
at J, ~ 1.45 for the two-dimensional KLM at half-filling 



24|. 



If the AF coupling between the local- 
ized spins is added, we favor antiferromagnetism even 
more. For this reason, J c is expected to increase with 
Jaf- Here, we do not attempt to provide the location 
of Jc(Jaf) for the two-dimensional HKLM, which is a 
formidable task, but only show the dominant tendencies 
in the problem. In Fig. 3(d), we present the intensity 
of the spin structure factor S(q) at <f=(7r, 7r) vs. Jaf, 
for several clusters at J=2.0 and half-filling Clearly from 
this figure, >!5(7r,7r) increases with Jaf, suggesting that 
J c will increase as well, as anticipated. As in the two-leg 
ladder [25| . we also have found no evidence of ferromag- 
netism close to half- filling for the 4x4 cluster. Note also 
that it is expected that the AF-phase will survive away 
from half-filling, for the three-dimensional KLM [3]. We 
believe this phase will also exist in the two-dimensional 
KLM close to half-filling. 

Discussion. Based on the numerical results described 
in this investigation, in Fig. 4(a) a qualitative phase di- 
agram for the HKLM close to half-filling is presented. 
The robustness of our results with respect to the shape 
of the cluster used, suggest our conclusions may be qual- 
itatively valid even in higher dimensions. The only mild 
assumption in Fig. 4 is the existence of an overlap be- 
tween the AF and d-wave superconducting regions. At a 
fixed J and with increasing Jaf this overlap is to be ex- 
pected. In the SC phase not overlapping with AF, strong 
short-range AF fluctuations are also to be expected. 

The phase diagram in Fig. 4 is in qualitative agreement 
with experimental results reported for some HF materi- 
als, such as Celn 3 and CePd2Si2, at low temperatures. 
A possible experimental trajectory is shown as a dashed 
line in Fig. 4. At ambient pressure these compounds 
are known to be AF [4]. As observed in the figure, for 
small values of J and Jaf (corresponding to low pres- 
sures in the experiments) there is no binding of holes 
(see Fig. 1) and superconductivity is not expected. How- 
ever, with increasing pressure at a critical value p Cl , and 
within a narrow pressure range, superconductivity devel- 
ops and coexists with long-range AF order This is 
compatible with our results, since at intermediate val- 
ues of J and Jaf (intermediate pressures) a tendency 
to d-wave superconductivity was numerically observed. 
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FIG. 4: (Color online) Schematic phase diagram of the HKLM 
close to half-filling. The solid line defines the region where d- 
wave SC should exist based on our numerical calculations. 
The yellow region is the long-range AF phase. The dashed 
line describes a trajectory with increasing pressure that is 
compatible with HF experiments (see text). 



At pressures even higher p C2 >p cl , the real HF system 
first stabilizes a SC state without AF long-range order, 
and then finally a transition to a non-SC paramagnetic 
phase is reached. Again, from the theory perspective 
this is reasonable since for even larger values of J and 
Jaf (larger pressures) the system eventually transitions 
into a phase without hole binding (see Figs. 1,2). Note 
that the positive slope of the proposed "path" in Fig. 4 is 
qualitatively correct since under pressure both J and Jaf 
increase, due to the increasing overlaps of wave functions. 
Note also that theoretically the maximum binding energy 
is reached at intermediate couplings in the dashed-linc 
path, suggesting that the SC critical temperature first 
increases under pressure, reaches a maximum, and then 
decreases, as in experiments ji, 26|. 

Summarizing, we have investigated the HKLM with 
emphasis on hole-pair formation and SC correlations. A 
region with a robust tendency toward a SC state was 
identified at intermediate values of J and Jaf close to 
half-filling. The binding of holes only appears when the 
AF interaction between the localized spins is considered, 
showing its critical importance for a proper theoretical 
description of heavy fermion materials. 
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